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Abstract 

A formulation is presented for modeling a resistive card in the 
context of the finite element method. The appropriate variational 
function is derived and for validation purposes results are presented 
for the scattering by a metal-backed cavity loaded with a resistive 
card. 
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1 Introduction 

A resistive card is an infinitesimally thin sheet of material which allows par- 
tial penetration of the electromagnetic field. Thin dielectric layers and very 
thin conductors whose thickness is less than the skin depth are examples of 
materials which can be modelled by resistive cards or sheets. Resistive cards 
are often used for radar cross section and RF power penetration control and 
as a result they have been studied extensively. Such studies have generally 
been done in the context of high frequency [1,2] and integral equation solu- 
tions [3,4], but to date the treatment of resistive cards within the context of 
the finite element method (FEM) has not been considered. Over the past 
few years, FEM has been applied to a variety of electromagnetic applications 
and it is thus important to incorporate the modeling of resistive cards in the 
FEM. In this paper we propose an FEM formulation which accounts for the 
presence of resistive sheets. To validate this formulation, results based on a 
physical modeling of the resistive sheet are also presented. In this case, the 
resistive sheet is equivalently replaced by a thin dielectric layer. The mod- 
eling of such a layer in the usual manner leads to larger and consequently 
inefficient linear systems, which is the primary reason for resorting to a math- 
ematical modeling of the resistive sheet. Results based on the mathematical 
and physical modeling are presented in connection with the scattering by a 
metal-backed cavity in a ground plane and these are used to validate the 
proposed mathematical model. 


2 Formulation 

A resistive sheet is characterized by its resistivity R which is measured in 
Ohms per square. Mathematically, it satisfies the boundary condition [5 

n x (n x E) = -Rh x (H + - H“) (1) 

where denotes the magnetic field above and below the sheet, E is the 
electric field and its tangential component is continuous across the sheet, and 
n denotes the unit normal to the sheet pointing in the upward direction ( + 
side). To a first order, this boundary condition can be used to simulate the 


2 


presence of a thin dielectric layer by setting [4,5] 


jk 0 (e r -l)t 

In this, t is the thickness of the layer, Z Q and fc 0 denote the free space intrin- 
sic impedance and wave number, respectively, and t T is the relative dielectric 
constant of the layer. Alternatively, a resistive sheet may be equivalently re- 
placed by a thin dielectric layer having thickness t and a relative permittivity 
of 


e, = I - 


jZq 

kjR 


( 2 ) 


Generally, the accuracy of this simulation increases as the thickness t is 
decreased. Typically, t should not exceed one-tenth of the wavelength in the 
material. 

Let us now consider a finite element solution of the fields within a volume 
V subject to a given excitation. The volume consists of some inhomogeneous 
dielectric having relative permittivity and permeability e r and fi T , respec- 
tively, and we shall also assume that resistive cards may be embedded within 
the dielectric (see Figure 1). In accordance with the finite element method, 
the volume is subdivided into M smaller volume elements and in this case 
we require that the resistive cards are tangential to the boundary surface 
of these elements. A weak solution of the fields within the volume can be 
obtained by extremizing the functional 


M 

f = £f', 


C= 1 


(3) 


F ' = J J J y P-(V x E) • (V x E) - k]t T E • E 

~\-jk 0 Z 0 #T'< H x h c )ds 


dv 


(4) 


with respect to the electric field E including that implied in H. In this 
expression, V e is the volume of the element which is enclosed by the surface 
S c and h c denotes the outward normal to S c . 
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Generally, for a dielectric volume not enclosing resistive sheets or other 
current sheets, the contributions of the surface integrals m (3)-(4) vanish 
everywhere except when S e coincides with the outer surface S a of the volume 
V This is a consequence of field continuity across the elements, but it a 
portion of the element’s surface coincides with a resistive sheet, then the 
surface integral in (4) does not vanish since the magnetic field is discontinuous 
as described in (1). Let us for example consider the surface S re which borders 
the eth and (e + l)th elements, and is coincident with a resistive sheet o 
resistivity R. Then the contribution from this surface to the surface integral 

in is r r t \ 

jk 0 Z 0 / / E- (H- x h rt )ds 
J J Srt 

from the eth element and 

-jk 0 Z 0 I / E • (H + x n re ) ds 

J J Sre 

from the (e + l)th element with n re pointing from the eth to the (e + l)th 
element. Combining these two integrals and employing (1), it follows that 
the contribution of the surface S re to the surface integral in (4) is 

- j k 0 Z 0 If [E • (H + - H ) x n re ] ds 

J J Srt 

= jk 0 Z 0 J J s J (nre X E) • {fl re X E) ds (5) 

Consequently, the functional F may be rewritten as 

f - 1 ///,R (V x E) • (V x E) - k 2 Q e T E • e| dv 

+jk 0 Z 0 J J s j (hr x E) • (h T x E) ds 

+jk 0 Z 0 (jj> s E- (H x h a )ds ( 6 ) 

in which S r denotes the surface occupied by the resistive sheet and S a is the 
outer surface enclosing the volume V. As usual, h T is the outward normal to 
S T and h 0 is correspondingly the outward normal to S a . If S T borders the 
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outer surface of the volume V, then S 0 should be considered to be just over 
the exterior side of S r (i.e., S 0 always encloses S r ). 

Having derived the explicit form of the functional F, we may now expand 
the element field using the standard linear shape functions. If the sources 
are within V then F should be modified to read 


F 



k 2 0 t r E E )dv 



+jk 0 Z 0 j j^(h r x E) • (h r x E )ds 
+j k 0 Z 0 ( j fj) E • (H x h 0 )ds 


( 7 ) 


where (J m ‘,M m ‘) denote the impressed sources internal to V. Then, upon 
setting the first order variation of F to zero, we can obtain a system of 
equations for the solution of the interior and boundary electric fields. For a 
unique solution of this system we must, however, specify a relation between 
the tangential electric and magnetic field which appear in the surface integral 
over S 0 - If we assume that the subject volume is that occupied by the metal- 
backed cavity recessed in a ground plane, as shown in figure 2, then S 0 reduces 
to the aperture area of that cavity. By invoking image theory, the magnetic 
field on the aperture can then be expressed as 


H = H* 


+Ir ''-W/ S „[( 




I+-VV]G 0 (r,r') 


[E(r) x z\ds' (8) 


where S a denotes the aperture surface, G 0 (r, r') is the free space scalar 
Green’s function, r specifies the observation point located on S a and I = 
xx + yy + zz is the unit dyad. Also, H ,nc denotes the magnetic field gen- 
erated by sources in the free space and H rc ^ is the corresponding reflected 
field when the cavity’s aperture is shorted. Substituting (8) into (7) gives 
a functional only in terms of the electric field. The system obtained from 
this functional will be partly sparse and partly full. In particular the volume 
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integrals and that over S r in (7) lead to a sparse submatrix involving the 
interior Helds of the cavity. However, in view of (8) the last surface integral 
of (7) over S c (or S a ) renders a full Toeplitz submatrix since the boundary 
integral is convolutional. Consequently, by resorting to an iterative solu- 
tion such as the conjugate or biconjugate gradient method in conjunction 
with the fast Fourier transform, the need to generate the Toeplitz matrix is 
eliminated thus maintaining the 0(n) storage requirement, characteristic of 
finite element solutions. The details pertaining to this implementation are 
discussed in [6] - [8]. In the next section we only present some results aimed 
at evaluating the accuracy of the proposed resistive sheet model. 

3 Numerical Results 

Let us consider the metal-backed rectangular cavity illustrated in Figure 2. 
The cavity is assumed to be empty (no internal sources) and is illuminated 
by a plane wave in the = 0 plane. For implementing the aforementioned 
solution, the cavity is subdivided into rectangular bricks and the results of the 
solution are shown in Figure 3. These are radar cross section (RCS) patterns 
and refer to a 1A deep cavity whose aperture is also 1A x 1A. The RCS 
pattern in Figure 3(b) applies to the cavity which is loaded with a resistive 
sheet of 100ff/D placed at its aperture, whereas the result in Figure 3(a) is for 
the untreated empty cavity. The simulation of the resistive sheet was done 
through direct discretization of the first order variation of the functional 
F as given in (7) and alternatively by modeling the resistive sheet as a 
dielectric layer of thickness A/20 having the dielectric constant computed 
from (2). As shown in Figure 3(b) the results based on the two simulations 
are in reasonable agreement and the differences among them is due to the 
finite thickness which was necessarily introduced in the physical model of the 
resistive sheet. As noted in [9] and [10], the dielectric layer introduces vertical 
components of the electric field which are not present in the resistive sheet. 
It is certainly of interest to point out that the presence of the resistive sheet 
at the aperture surface reduced the RCS of the cavity by lOdB at normal 
incidence and by as much as 20dB at grazing incidence. 

The second geometry which was considered is a circular metal- backed 
cavity again situated in a ground plane. In this case the aperture of the 
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circular cavity is loaded with a sheet having non-uniform resistivity given by 


R(p) = 


0.1 + 10 (^) : 


t OO 


p > 0.5a 
elsewhere on S a 


( 9 ) 


where a denotes the radius of the aperture. Results with and without resistive 
loading for a cavity having a = lin. and a depth of 0.25in. are shown in 
Figure 4. These RCS patterns were computed at 16GHz and the incident 
field was a plane wave polarized along the 6 or 4> directions. Again, the data 
in Figure 4 demonstrate the validity of the proposed mathematical model. 
Also, as in the case of the rectangular cavity the presence of the resistive 
cards reduces the overall RCS of the cavity and this reduction is primarily 
due to the recced field intensity near its perimeter. 


4 Conclusions 

A formulation was derived for modeling resistive cards within the context of 
the finite element method. Essentially, the pertinent variational functional 
was supplemented with an additional boundary integral over the surface of 
the resistive sheet/card. Results based on the discretization of the functional 
were also presented and these were aimed at demonstrating the accuracy of 
the derived mathematical model. 


References 

[1] T.B.A. Senior, “Scattering by Resistive Strips,” Radio Science, Vol. 14, 
pp. 911-924, 1979. 

[2] M.I. Herman and J.L. Volakis, “High Frequency Scattering by a Resis- 
tive Strip and Extensions to Conductive and Impedance Strips,” Radio 
Science, Vol. 22, pp. 335-349, May-June 1987. 

[3] E.F. Knott and T.B.A. Senior, “Non-Specular Radar Cross Section 

Study,” T niversity of Michigan Radiation Laboratory Technical Report 
Oil 764- x (also U.S. Air Force Avionics Lab. report AFAL-TR-73-422) 

January 1974. 


7 



[4] R.F. Harrington and J.R. Mautz, “An Impedance Sheet Approximation 

for Thin Dielectric Shells,” IEEE Trans. Antennas Propagat ., Vol. AP- 
23, pp. 531-534, 1975. 

[5] T.B.A. Senior, “Combined Resistive and Conductive Sheets, IEEE 
Trans. Antennas Propagat., Vol. AP-33, pp. 577-579, 1985. 

[6] J.M. Jin and J.L. Volakis, “A Finite Element- Boundary Integral Formu- 
lation for Scattering by Three-Dimensional Cavity-Backed Apertures,” 
IEEE Trans. Antennas Propagat., Vol. AP-39, pp. 97-104, January 
1991. 

[7] J.M. Jin and J.L. Volakis, “Electromagnetic Scattering by and Trans- 
mission Through a Three-Dimensional Slot in a Thick Conducting 
Plane,” IEEE Trans. Antennas Propagat., Vol. AP-39, pp. 543-550, 
April 1991. 

[8] J.M. Jin and J.L. Volakis, “Scattering and Radiation From Microstrip 
Patch Antennas and Arrays Residing in a Cavity,” submitted to IEEE 
Trans. Antennas Propagat.. 

[9] T.B.A. Senior and J.L. Volakis, “Sheet Simulation of Dielectric Layers, 

Radio Science , 22, Dec. 1987, pp. 1261-1272. 

[10] J.L. Volakis, “Numerical Implementation of Generalized Impedance 
Boundary Conditions,” 1989 URSI Electromagnetic Theory Symposium, 
Aug. 1989, Stockholm, Sweden. Symposium Digest pp. 434-437. 


8 



FIGURE CAPTIONS 


Fig. 1 Cross Section of a dielectric volume enclosing a resistive card (a) card within the 
dielectric (b) card on the surface of the dielectric 


Fig. 2 Geometry of a cavity-backed aperture in a ground plane. 


Fig. 3 Monostatic RCS for a rectangular aperture (R x R) backed by a rectangular cavity 

(R x R x R) in the <t>=0 plane, (a) Empty cavity without resistive loading (b) 
Empty cavity whose aperture is covered with a resistive sheet having a resistivity of 

100Q/n. (Solid and dashed lines correspond to results based on the mathematical 
simulation; circles and squares refer to results based on the physical modeling of 
the resistive sheet) 


Fig. 4 Monostatic RCS for a circular aperture (1 inch in diameter) backed by a circular 

cavity (1 inch in diameter and 0.25 inches deep) at 16 GHz. (a) Empty cavity, (b) 
Empty cavity covered with a resistive sheet whose resistivity is given by equ. (a). 
(Solid and dashed lines correspond to results based on the mathematical 
simulation; circles and squares refer to results based on the physical modeling of 
the resistive sheet) 




Figure 2 
















